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Abstract 

In a recent paper we have introduced several possible inequivalent descriptions of the 
dynamics and of the transition probabilities of a quantum system when its Hamiltonian 
is not self-adjoint. Our analysis was carried out in finite dimensional Hilbert spaces. 
This is useful, but quite restrictive since many physically relevant quantum systems live 
in infinite dimensional Hilbert spaces. In this paper we consider this situation, and we 
discuss some applications to well known models, introduced in the literature in recent 
years: the extended harmonic oscillator, the Swanson model and a generalized version of 
the Landau levels Hamiltonian. Not surprisingly we will find new interesting features not 
previously found in finite dimensional Hilbert spaces, useful for a deeper comprehension 
of this kind of physical systems. 
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1 Introduction 


In ordinary quantum mechanics one of the fundamental axiom of the whole theory is that the 
Hamiltonian H of the physical system is self-adjoint: H = W. This condition, shared also 
by all the observables of the system, is important since it ensures that the eigenvalues of 
these observables, and of the Hamiltonian in particular, are real quantities. Moreover, since 
the time evolution deduced out of H is unitary, it preserves the total probability: if T(t) is a 
solution of the Schrodinger equation then ||T(t)|p does not depend on time. 

This is clear since, if H does not depend explicitly on time, and is 

unitary, hence norm-preserving. Of course, this is false ii H ^ and in fact, in this case, 
||'h(f)||^ does indeed depend on time, in general. Sometimes this is exactly what one looks 
for: in many simple systems in quantum optics, for instance, non self-adjoint Hamiltonians 
are used to describe some decay, so that there is no reason for the probability to be preserved 
in time. In other situations, one would prefer to avoid any damping, so that the aim is to 
hnd some way to recover unitarity even when H ^ and in fact several attempts have been 
proposed along the years by different authors to discuss this and other aspects of time evolution 
for systems driven by non self-adjoint Hamiltonians. Here we refer to n-i, and references 
therein. We also suggest PCQ] for two rather general, but not so recent, reviews on this and 
related subjects, and for a more recent volume, rather mathematically oriented. Recently, 
in [12] , we have discussed some dynamical aspects of this kind of systems, working always with 
hnite-dimensional Hilbert spaces, to avoid dealing with unbounded operators and to use the 
elegant mathematics of pseudo-fermions (PFs), [H] [HI [T5| , in the analysis of the systems. The 
aim of this paper is to discuss what should be changed and considered with more attention 
when going from a finite to an infinite dimensional Hilbert space. We will see that the general 
structure survives to this transition, but that, not surprisingly, several mathematical subtleties 
must be properly taken into account. Also, we will discuss how this transition produces non 
trivial physical consequences. 

The paper is organized as follows: 

In the next section we review the general functional structure associated to a non self-adjoint 
Hamiltonian, and its dynamics, as well as few inequivalent dehnitions of transition probability 
functions, recalling what was done in [12]. In particular, we discuss what can happens in 
presence of a metric operator. In Section [3] we apply our general results to some well known 
(and reasonably simple) models, the extended quantum hamonic oscillator, the Swanson model, 
and a generalized version of the Landau levels, see [16] for a recent review on these models. 
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Our conclusion is that it may be not so convenient to adopt a different scalar product, 
the one which makes oi H a. self-adjoint operator, even if this is possible, in many explicit 
situations. The reason is that, as we will see in Sections l3.2l and [373l this will force us to somehow 
change the original model, by imposing extra constraints on the coefficients appearing in the 
dehnition of the model, coefficients which, in our examples, measure the non self-adjointness of 
the Hamiltonian. 


2 The general settings for H ^ W 


As we have already said, in this paper we will mainly be interested in inhnite dimensional 
Hilbert spaces. Then, our operators can be unbounded, as it very often happens in quantum 
mechanics, even for very simple system^. The main ingredient is an operator H, acting on 
the Hilbert space T-L, with H ^ H\ and with all (multiplicity one) real distinct eigenvalues 
En, n = 0,1,2,.... This condition is useful only to simplify the notation, and could easily 
be revomed. Here, the adjoint W of H is the usual one, i.e. the one dehned in terms of the 
natural scalar product (.,.) of the Hilbert space "H: {Xf,g) = (/, for all f,g E Ti such 

that / G D{X) and g G Zl(Xl), the domains of X and Afl respectively. We call ||.|| the norm 
dehned by (.,.). 

In this paper we will consider also two other scalar products on "H, (., .)^ and (., The 
reason for that will be clarihed in the following but can be simply understood already at this 
level; even if {Hf,g) ^ {f,Hg) in general, it will happen that {Hf,g)^ = {f,Hg)^ and that 
9) = \f 1 H^9 / • Of course, due to the fact that H may be unbounded, we have to take 


f,gE D{H) and f,gE D{W), and these two sets must be, if H and W are unbounded, dense 
subsets of H. Needless to say, nothing like this occurs for hnite dimensional Hilbert spaces, 
which is the case most of the time considered in the literature, see for instance mm- in fact, 
in this case, both D{H) and D{H^) coincide with the whole Hilbert space. 

We consider here an Hamiltonian H with distinct real eigenvalues, corresponding to different 
eigenvectors (p*, G "H, A; = 0,1, 2, 3,...: 


H^pk E]«Lpi~. ( 2 . 1 ) 

The set = {<pk, A: = 0,1,2, 3...}, in many papers, is assumed to be a basis for "H. This is 

^It is well known, in fact, that the Hamiltonian H of one of the simplest, and more important, quantum 
mechanical system, the harmonic oscillator, is unbounded, as well as the lowering and raising operators in terms 
of which H can be factorized. 


3 








surely true if dim(?^) < oo, because the various vectors are automatically linearly independent, 
since they correspond to different eigenvalues. However, in several recent (and not so recent) 
studies, this has been shown not to be true in general, see [16] and references therein, when 
dim('H) = oo. Therefore the assumption that is a basis for "H, if dim{'H) = oo, is not 
entirely justihed. On the other hand, in all the systems in our knowledge, what is true is that 
is a complete set in "H. This is, of course, a big mathematical differenci^, [H] . For this 
reason, we will restrict here to this lighter, and always satished, assumption; is complete in 
'H. To this working assumptions we have to add similar conditions concerning iff which, as an 
operator on "H, may appear very different from H. In particular we will assume here that 
admits eigenvectors with the same eigenvalues as H: 

= EfcTfc. (2.2) 

Hence H and W are assumed to be isospectrals, in this paper. The set = {Tfc, A; = 
0,1, 2, 3 ...} will also be assumed to be complete in "H, but not necessarily a basis. Of course, 
together with these are biorthogonal sets: ((p„, = 6n^k- 

Similarly to what is proposed in [T8|, rather than working on "H with the complete sets 
and it is convenient to dehne a subset of "H, which we call Vp^ys, as follows: 


V. 


phys 


Ifen-. f = \ 

\ n n ) 


(2,3) 


Then, Vyhyg is not necessarily all of H, except, for instance, if l-L is finite dimensional or if 
and are biorthogonal (Riesz) bases. The subscript phys stands for physical, meaning with 
this that the set Vphys is assumed to contain all the physically relevant vectors needed for the 
full description of the physical system S we are interested in. In other words, even if Ty, and 
are not bases for H, they are still sufficient to expand all those vectors which have a physical 
meaning. Of course, an obvious requirement is that each and each \l/„ belong to Vphys, since 
these vectors are surely the most physical of the system, being eigenstates of H and i/b Due 
to the biorthogonality of Vp and this is equivalent to require that, for all k, 

(<dn, V?fc) (2.4) 


set X = {xn, n G N} is complete in H if, for any f G H and for any e > 0, it exists > 0 and 
a set of (complex) coefficients k = 1,2,..., N, such that ||/ — < £■ In general these 

coefficients depend on N. When it happens that they do not depend on N, the set X is a basis for TL since 
11/ — CkXkW —>■ 0 for N diverging. It is well known that completeness of X is equivalent to X being a 

basis if the vectors x^s are mutually orthonormal. Otherwise this is not true: any basis (orthonormal or not) 
is complete, but not all complete sets are bases. 
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Of course, these equalities surely hold whenever and collapse to a single orthonormal 
(o.n.) basis of "H, when they are Riesz bases, or for hnite-dimensional Hilbert spaces. In other 
cases, however, they must be explicitly checked. Notice that, since and are complete in 
"^5 T^'phys is dense in "H. Then, Vp^ys is a large set, indeed. 

Remark:— It could be interesting to notice that, while here the physical space is dehned 
by Ty, and which in turns are determined by H and Ri, in [18] the Hilbert space where 
the model lives has to do with the set of all the observables of the system. This is a natural 
procedure, of course, even if it is not guaranteed that the intersection of all these domains is 
a sufficiently large set, in concrete situations. We also would like to notice that, contrarily to 
what stated in [18], the analysis carried out in this paper will suggest that the natural scalar 
product to be considered in the physical space is exactly the one originally dehned in "H. This 
will be clarihed in the second part of the paper. 

In [12] we have discussed how the dynamics of a physical system S should be dehned when 
its Hamiltonian H is not self-adjoint. Our suggestion, which agrees with the point of view 
widely adopted in the literature, is that the wave function <l)(t) of S should satisfy a standard 
Schrodinger equation, i^{t) = with $(0) = $ 0 . Of course, it is natural to require hrst 

that $0 ^ 'I^phys, and to assume also that this property is preserved under time evolution. In 
other words, we would like to have <h(t) G Vphys for all t > 0, and not just for f = 0. This is not 
granted, because it is not necessarily true that, even if $0 = Xln *ho) *^ 0 ) 

then 

m = = E (2.5) 

n n 

However, this is the case if, for instance, (z) is bounded and (m) ((p„,<h(t)) is a 

Cauchy sequenc^ in "H. In fact, under these assumptions we have: 

$(t) = (T„, $ 0 ) (p. 

V n 

using the continuity of Hence 

® 0 > = d>o) = 

n n n 

n n 

^Of course, this is surely true for t = 0. We are here assuming that this is also true for t > 0. 
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which is half of equation (I2.5p . As for the second half, since ^(t)) is, by assumption, 

a Cauchy sequence for all t > 0, it surely converges to some vector of l-i. The fact that this 
vector is exactly $(f) follows from the fact that is complete in "H and by the biorthogonality 
of and since we have ((p^, [d’(t) — 'Ylm (V’n, ^(t)) ^n]) = 0 for all k. 

This result shows that, when the dynamics of the physical system S is driven by a non 
self-adjoint Hamiltonian H, problems arise both at a pure algebraic level and at the level of its 
dynamical description. If, from one side, the introduction of Vphys looks quite reasonable on 
a physical ground, requiring that Vp^ys is stable under time evolution is not granted a priori, 
and some extra conditions are required. Of course, the ones given here are sufficient conditions, 
so we can imagine that they can be lightened. Once again, all the conditions are satished if 
dim('H) < oo, which is the situation quite often discussed in the literature, or when and 
are biorthogonal Riesz bases. 

Remarks:— (1) Of course, these steps can be slightly modihed, with few and obvious 
changes, if one assumes as driving Hamiltonian W rather than H. 

(2) The dynamics of the operators in the Heisenberg, as well as in the Schrodinger, repre¬ 
sentation is not uniquely defined. A natural choice is 

X{t) = (2.6) 

As we have discussed in [12], this is not the only possibility, and it is not necessarily the 
most convenient, since neither nor are unitary (which however is exactly what one 
looks for, sometimes). Moreover, adopting the rule in (12. 6 h it is not so easy to find integrals 
of motion for the system, since [H,X] = 0 does not imply that X{t) = ^(0) for all t > 0. 
Finally, a serious difficulty is that the time evolution is no longer an automorphism of the set 
of observables, since in general {XY){t) ^ X{t)Y{t), and this complicates in an enormous way 
all the computations. We refer to na for further considerations on this aspect, which is not 
our main concern here. 

(3) It might be more convenient to replace the set Vphys with 

'^phys-=lf^'^-f = ^{Vn,f)'^n, Or / = ^ (^'„,/) (/? 

n n 

Of course, each element of Vphys also belongs to while the vice-versa is not true. This 

means, in particular, that is dense in T-L. We prefer Vphys because of a more evident 

symmetry between Xy, and X^,. 
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As in na, we are mainly interested here in defining the probability transition between two 
states, the initial state of the physical system S, $ 0 , and the final state, $/. As we have 
discussed in [12], this definition is not unique. This is due to the presence of, at least, three 
inequivalent scalar products defined in Vp^ys or even in "H. In analogy with [12], we call these 
products and and we call b and H the related adjoint maps: {Xf,g)^ = 

and {Xf,g)^ = {f,X'’g)^, for all / and g in H and for all operators X for which 
these equalities make sense. The products (., .),j, and (., .)^ are such that, as discussed before, 
and H are self-adjoint with respect to them: 

(/, Hg)^ = {Hf, g)^ , (/ft/, ^ = (/, H^~g^^ , (2.7) 

for all f,g E D{H) and f,g& D{W). Here we are implicitly assuming that both D{H) and 
D(i/t) are subsets of Vphys- Using (I2.7jl we deduce that H = H\ and that i/t = (f/tj^. Once 
we have three scalar products, we also have several different possible definitions of the transition 
probabilities. The ones we consider here are the following: 

(■K/. ■!>(*)>,. ' 

( 2 . 8 ) 

These were already introduced and analyzed in [12], for systems living in finite dimensional 
Hilbert spaces, and we have seen that they produce different results, so that, in fact, they are not 
physically equivalent at all. Then we have proposed to discriminate among these three defini¬ 
tions using some concrete experiment, and we have considered a simple two-level system. Here 
we want to carry on a similar analysis, but considering systems which live in an infinite di¬ 
mensional Hilbert space. This will be done in Section [3l where we will see how it is possible, 
in principle, to discriminate among the functions in fl2.8p . in order to understand which is the 
most appropriate expression of the transition probability, and why. Incidentally we observe 
that, with these definitions, the images of and is always the 

set [0,1], for all f > 0. 

Remark:— If we were interested in keeping the time evolution unitary, then would 

be the more natural choice, since H is self-adjoint with respect to (., .)^, and therefore 
is unitary with respect to this scalar product. However, here we are much more interested in a 
comparison between the theoretical results with some experimental data. This is also relevant 
in view of the fact that, if we replace H with the, equally valid, operator H\ the unitarity 


P^Q —* 








(t) := 






P‘p 


(f): = 
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requirement would suggest, of course, to use In other words: different ingredients 

produce different rules. 

Interestingly enough, we will see that the analysis of the functions in (I2.8j) in some concrete 
models suggests to avoid the use of (., .)^ and and to restrict to and to 

as a consequence. The other choices, in fact, sometimes produce a unwanted extra constraint 
in the range of the parameters of the model, in order to make sense out of the model itself. 
This will be made explicit in Sections 13.21 and 13.31 


2.1 Refining the structure 

What we have discussed so far does not imply the existence of any particular relation between 
the three scalar products above. They are just related, in principle, to H and HK This is 
because the two new scalar products are introduced here just to make of H and two self- 
adjoint operators. In many physical systems considered in the literature, however, a relation 
between them does in fact exist, and it is provided by the so-called metric operator. Again, 
while there is no problem to introduce this operator if dim('H) < oo, serious problems may occur 
for inhnite-dimensional Hilbert spaces. The reason is that it may happen that this operator, or 
its inverse, or both, are unbounded. When this happens, we have to pay attention to domains. 
In particular, we dehne 


D{S^) := 


|/ G P : {ifn, f) G , 


D{S^) ;= g 


GP; 

n J 


and 


^ 1 /) ^n- P’^9 ni g) ^ ni 


(2.9) 


for all / G P(5'^) and g G D{Sqi). To make the situation technically simpler, it is convenient 

to work under the assumption that Vphys ^ fl D{S^). This makes of and two 

densely dehned operators, if Vphys is dense in P, as we have observed in several concrete 

examples discussed so far, [16]. These operators have the following properties: (i) G D{Sp), 

and = (pn] (i) Pn G D{S^,), and (iii) they are positive operators, and, under 

suitable conditions, they admit a self-adjoint (Friedrichs) extension, which are also positive, 

1 /2 

and which we indicate with the same symbols; (iv) these extensions admit square roots, Sp 
1 /2 

and ; (v) for all f,gE Vphys the scalar products introduced above are related as follows: 


(/,9)i. = (/.S»9> = (sF/.SF9), 






and 


U,9)^={f,S^a) = {Sy^f,Sy^g). 

Another feature of S^p and is that, again under suitable assumptions, they relate the different 
adjoint introduced so far, f, b and H. In fact, taken an operator X of iS, and assuming for 
simplicity that X leaves invariant 'Dphyg together with X\ X^ and X**, we deduce the following 
equalities: 

X^f = S^X^SJ, X»/ = S^X^S^f and X^/ = SlX^S^f, 

for all / G Vphyg. Of course, we are also assuming that Sy, and 5*^ leave Vp^yg invariant. Finally 
a direct computation shows that 

Remark:— The results sketched in this section suggests that it is the Hamiltonian H itself 
which somehow hxes its preferred Hilbert space. This is because both (., .)^ and (., .)^ are 
dehned via Sp and which are constructed, in turns, by the eigenvectors of H and W. This 
is similar to what happens in algebraic quantum dynamics, see [I9] and references therein, 
where the Hamiltonian (self-adjoint, in that context) is used to dehne a suitable topology on 
the algebra of the operators needed in the description of the physical system. This aspect is 
also discussed in na, together with the role of non zero temperature states. 

If we now introduce formally Hq := and we see that 

Hq = hI = W, and that 8 = {e^, k = 0,1,2,...,} is an orthonormal (o.n.) set of 71 
of eigenstates of Hq: HoCk = EkCk- To move from formal to rigorous results we need to perform, 
of course, a deeper analysis of the operators considered. In particular, if for some reason Vphys 
is left invariant by H, and >5^^, as it happens in some concrete examples, [16], then Hq 
turns out to be a densely dehned symmetric operator and, if Jy, is a Riesz basis, then 8 is an 
o.n. basis. 

Remark:— A possible alternative approach consists in introducing, together with Vphys, the 
set of what we can call physically relevant operators, Ophys, as the set of all the operators X, 
bounded or not, densely dehned on "H, which leave stable Vphys together with all their adjoints, 
Xt, X** and X^ Then, our working assumption is that H, and belong to Ophys- 

This is not very diherent from what it is done in the literature on unbounded operator algebras 
when one introduces the set Cl{V), which is the *-algebra of all the closable operators dehned 
on the dense set V which, together with their adjoints, map V into itself, [19]. In fact, also in 
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our case, the set Ophys turns out to be an algebra of unbounded operators, having Vphys as the 
common domain. One of the obvious differences between Ophys and 0'{T>) is that the hrst one 
involves several involutions, while the latter just on^. 

3 Examples 

In this section we will consider some concrete examples, already considered in the literature, to 
compare the different dehnitions of transition probabilities introduced in (I2.8jl and to deduce in 
this way which one, among the different possibilities, is the more appropriate. More concretely, 
we will analyze what happens for the extended quantum harmonic oscillator (EQHO), for the 
Swanson model, and for the extended Landau levels (ELLs), see [20l [211 122], computing the 
different expressions of the transition probabilities for particular choices of the initial and the 
hnal states. With respect to what has been discussed in na, the role of unbounded operators 
will clearly show up, and we will see that new phenomena will occur exactly because of the 
inhnite dimensionality of the Hilbert space, leading to conclusions somewhat different from 
those deduced in [T2] . 

3.1 The extended quantum harmonic oscillator 

The Hamiltonian of the EQHO, as originally proposed in [20] and then rewritten in terms of 
pseudo-bosonic operators in [21], looks as follows: 

^ + *\/2p, 

where z/ is a strictly positive parameter and [x,p] = -il. Here II is the identity operator on 
Pi = 12^ (M). Of course, is manifestly non hermitian. 

Introducing the standard bosonic operators o = (x -|- ^), — ^), [a, al] = 1, 

and the related operators Ay = a — and By = aP + we can write Hy = v {ByAy yj, 1), 
where '^y = It is clear that, for aA v > Q, Al ^ By and [Ay, By] = 1. Hence we have to 

do, at least formally, with pseudo-bosonic operators. We refer to [1611201121] for more details. 
In particular, we have proven that the sets = {^pn\x)} and = {\l'i^^(x)} of eigenstates 
of Hy and Hi are not biorthogonal bases, but still they are "D-quasi base^, and that they are 
both complete in £^(M). 

"‘This is not really so, because one has also to deal with the restriction of the involution to V, but this has 
nothing to do what we are considering in this paper. 

®This means that and still resolve the identity, but only weakly on the dense set T). 
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Also, we have deduced that, with a proper choice of normalization. 




\ dx^ u ^ 


and 



( d 

_I ^ — _ 

yj2^ n\ Y dx 



g-i(x+V2/y)2^ 


They both correspond to the same eigenvalue, En'^ = i/{n + 7 ^), for and Hi respectively. 
The so-called metric operator 0,^, mapping into is a simple multiplication operator, 

which looks like 0 ,^ = e-^(“+A) _ gS/i/^ g- 2 ^/ 2 ^_ have T^r^(a:) = 0 y(^n^^(x), for all n. 

It is clear that 0y is unbounded, since it is not everywhere dehned in £^(M). However, it 
is invertible with unbounded inverse and we obviously have that {x) G D{Q^), and that 
G D{Q~^) for all n. Of course, the operator Oj, must be identihed with the operators 
5',^ introduced in fl2.9p . 

In order to compare the transition probabilities in fl2.8p . and to make the computations 
simpler, we restrict to the hrst two eigenstates of and Hi: 






A’U) = —Ah) = — 


and 


TflA 

7/A 


71 - 1/4 

1 /A 


yj-i/4 ^i/4 


It is now an easy exercise to compute the transition probabilities for some different choices of 
$0 and $/. For instance, if $0 = and $/ = we hnd that 


which are all different and independent on time. Notice that when u is taken large enough the 
three probabilities all converge to one. This appears in agreement with the fact that, in this 
limit, the self-adjoint part of Hi^ is much larger than the remaining part, so that, essentially, 
the deviation form a standard situation is really small. 

If we now take = $/ = it is easy to see that = 1, 

which is not surprising. In fact, this result can be easily generalized to the following situation: 
suppose that iS is a physical system with (non self-adjoint) Hamiltonian H and let pa be an 
eigenstate of H. If we choose $0 = <F/ = then = 1. 


11 














In this case the three probabilities coincide. Hence, this choice is not useful to discriminate 
among the various functions in (I2.8p . 

More interesting for us is the situation in which $o = and In this 

case, obviously, we have <I>(f) = and after some computations 

we conclude that ^ 

(t) 2 (2 + + 2 ncos(nf)) ’ 

((^ “ ^) cos{iyt))^ , 

and 

(P _ (1 -!)' + !(! + cosjut)) 

2 8 + i ^2 + 4j,cos(nf) ' 

These formulas show that each probability transition go to | when z/ —)■ oo, for all t. This is in 
agreement with our previous interpretation of this limit. Also, if /v —)■ 0 then all these functions 
converge to zero: the <I)o —)■ $/ transition is not allowed, in this case, whichever choice we do. 
Notice also that the three functions are periodic, with a period which is exactly the smaller 
the value of u, the longer the period. 

It is not difficult now to imagine, at least in principle, concrete experiments capable to 
discriminate among the three dehnitions in ( 12 . 8 p . just comparing the experimental results with 
what we have deduced above with the hrst or the third choice of <I>o and $/. The second choice 
is not useful for us, since the resulting functions do coincide. We notice that, in the analysis of 
the EQHO, there is no reason to prefer to the other two possibilities, expect for its 

agreement with experiments. We will see that this is not what happens for the other models 
we are going to consider next. 

3.2 The Swanson model 

It is interesting to discuss also what happens for the Swanson model, since, as we will see in a 
moment, these are new facts which were not observed in [ 12 ] and in the previous example, the 
EQHO. Again we refer to [I 6 ] for the details of our construction. 

The non self-adjoint Hamiltonian of the model is 

^0 = ^ (/ + ^ tan( 20 ) - x^) , 

where 6 * is a real parameter taking value in (—f, f) \ { 0 } =: /, [20 ], and [x,p] = il. 
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Introducing now the (standard bosonic) annihilation and creation operators a, a\ and their 
linear combinations 


Aq = cos{6) a + isin( 6 ') + e , 

Be = cos(d) i sin( 6 ') o = :^ , 


we can write He = (Be Ae + , where ue = dehned because cos(26*) 7 ^ 0 for 

all 6 E I. The eigenfunctions of Hq and hfj, forming the sets and have been found 
in |21j : 

ip^n\x) = Hn (e^^x) exp {-i x^} , 

{“i , 

where Hn{x) is the n-th Hermite polynomial. These functions all belong to £^(M), but they 
are not bases of this space, thought being complete, 


The operator S',!, which maps Jy 


(8) 


(O') 

into is the following: 


(Svt/) (x) = e (e , 


which is, of course, not everywhere dehned. Hence, is unbounded, with unbounded inverse 
satisfying the following: (S^f) (x) = (e^^x). It is easy to check that ipn\x) E ^^(S^r) 

and that S^v9 n^(x) = for all n > 0. Analogously, we can check that \l/n^(x) G D{S^) 

and that S,pTn^(x) = ipn\x). 

In the attempt to deduce some transition probability function explicitly dependent on time 
we consider hrst the initial state of S to be <I)o = (pQ +ip\ , and $/ = Tq , as in Section 
13.11 However, this is not enough, due to the fact that the vectors in and are just 
rotated versions of the eigenstates of a quantum harmonic oscillator, with a slightly different 
normalization. In fact we get 

^ cos^(2g) 

1 + cos2(20)’ 

which does not depend explicitly on time. A similar result can be deduced also for the other 
functions in fl 2 . 8 l) : 


pvt pv _ Vcos(4g) ^ (cos(4g))^/^ 

) 2 cos 26* ’ ^ cos26'(cos(40) + !)■ 

It is important to notice that these two last functions are not dehned for all 9 E I, but only 
if 6 * G (— f 1 f) \ {0} =: h- This suggests that the dehnitions and are 

somehow artihcial, and only make sense if we are willing to change the original model by 
imposing, as we have to do here, more constraints on the parameters of the system. 
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Remark:— It is interesting to observe that, when 0 —)■ 0, all the probabilities converge to 
as it is expected to happen: in fact, in this case. He becomes the Hamiltonian of the standard 
quantum harmonic oscillator. 


m 


The same conclusion can be deduced if we compute the three probability transition functions 

(Q\ /n\ /g\ /g\ 

taking $0 = = ^0 + • The main difference, in this case, is that 


an explicit dependence on time appears. In fact we get: 

2cos^(26*)(l + cos{ue)t) 


H 




{t) = 


(1 + cos(26'))^ 


and 

^ cos^/2(46')(l + cos 2(26') + 2cos(26') cos(a; 0 t)) 

- 2(1 + cos(40)) cos3(20) 

Except for which is dehned for 9 E I, and again make sense 

only if 6 ^ G Ii. Then the conclusion is the following: when S is driven by an Hamiltonian 
H which is not self-adjoint (and has real eigenvalues) it seems more convenient not to change 
scalar product, looking for some other scalar products which makes of P a self-adjoint operator. 
This is because, if we do that, we may need to restrict the original range of the parameters of 
the model modifying, in fact, the original model. 


3.3 Extended Landau levels 

A similar conclusion can be deduced by considering this third model, whose main ingredients 
are the operators dehned in [22] as follows: 

^ ^ (4 - tdy + f (1 + 2h) - f (1 - 2h)) , 

Bl = ^ ( — (9a; — idy -\- |(1 — 2 /^ 2 ) + ^(1 + 2/Cl)) , 

^2 = -^ ( —*(9a; -\- dy — y(1 -|- 2 /C 2 ) -|- |(1 — 2kl)^ , 

.^^ = 72 -dy + W- 2 ^ 2 ) + 1(1 + 2 ki)) , 

where kj G ]—They satisfy the two-dimensional pseudo-bosonic commutation rules 
[Aj, Bk] = The vacua of Aj and Pj are 

j (po,o(x,y) = Ny, exp |-x(l + 2 ^ 2 ) - - 2ki)^ 

\ vl>o,o(x, y) = exp {-^(1 - 2 ^: 2 ) - ^(1 + 2k ,)} , 

where Ny, and are normalization constants which are chosen in such a way that ( 9 ?o, 0 ) To,o) = 
1. We hx them as Ny, = Of course, since k, and k 2 are such that — 1 < /cj < |, 

j = 1,2, both these functions are square integrable. 
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As in 1221 we define the fnnctions 



where n,l = 0,1, 2, 3,..., and the biorthogonal sets n,i{x, y), n,l > 0} and = 

y), n,l > 0}. Introdncing fnrther hi = BiAi — and ^2 = B 2 A 2 — |l, it is clear that 

[hi, h 2 ] = 0, and that 



and 



Fnrthermore, dehning S^p = e and S,;, = can check that, for in¬ 


stance, Sq,ipn,i = '^n,i, for all n and 1 : in this case, the metric operator is a simple multiplication 
operator. Also, it is clearly unbounded with unbounded inverse. 

Remark:— The particular case ki = k 2 = 0 returns the standard Landau levels, [22] . This 
can be understood already from what we have discussed here: in this case, in fact, F'lp = S',f = 1, 


and the sets Jfo and collapse to a single set of o.n. functions, complete in 


Let now assume that the system evolves according to a very simple Hamiltonian: H = 
hi + = BiAi, and that it is prepared in the state $0 = v?o,o + 9 ^ 1,0 + ¥^ 0,1 • We want 

to compute the transition probabilities, as introduced in (12.Sh . to hnd S in the hnal state 
= d'0,0 + 

After some lengthy computations we hnd that 


= 2(1 cos{t))pki,k2 


where 


Pk\,k 2 • 


y/(l - 4fc2)3(l - 4fc2)3 


while the expression for 




(2 -|- 3ki — 3 / c 2 — 4:kik2){3 + Aki — Ak2 — ^kik2) ’ 

(t) looks more complicated, since it is not possible to separate 


the dependence on time and on kj in the hnal formula: 



(1 + 2A:i)(l -2A;2) 


1 + ^1 — k2 


^ -|-2cos(f) 


(l + 2A:i)(l - 2 ^: 2 ) j 


l + ki- k2 
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where 


_ (l + 4fei)(l-4fc2) _ 

9k^,k2 ■ ^ 2 A;i)(1 - 2k2){l + Ski - Sk2 - 8kik2)' 

A similar formula could also be deduced for Notice that we get = 

= ^(l+cos(f)) if we go back to the ordinary Landau levels, i.e. if we take ki = k 2 = 0. 
It should be emphasized here that (t) can be found only under the additional requirement 

that kj takes value only in ] — | [ rather than in the original, larger set. This is very close to 

what we have seen for the Swanson model, and again the suggestion is that, if we don’t like 
this kind of additional restrictions on the parameters of the model, we have to choose, among 
the possibilities given in (I2.8p . the original one, i.e. working with the scalar product 

(.,.): this is something new with respect with what was found in [12], and with what is quite 
often discussed in the literature, since is a phenomenon which can only be seen in an inhnite 
dimensional Hilbert space. 

4 Conclusions 

After a general discussion which extends to an inhnite dimensional framework what originally 
proposed for the dynamical problem generated by a non self-adjoint Hamiltonian acting on 
hnite dimensional Hilbert space, we have deduced some consequences of our choices in the 
computation of several, inequivalent, transition probabilities. With the help of three simple 
examples we have seen that, in order to keep unchanged the original features of the model 
under analysis, and in particular the ranges of the parameters dehning the model, the only 
possible choice of the transition probability is This, in turn, suggests that the only 

realistic scalar product is the original one, (.,.), while (., .),j and (., .)^ should be understood 
only as auxiliary useful tools in the analysis of the model, but not really essential or having 
any deep physical interpretation. 

Of course, this is not really so if we admit the possibility of changing the model on the way, 
i.e. to further restrict the values allowed for the parameters to ensure, as in this paper, the 
square-integrability of some relevant functions to respect to some particular metric. In this 
case, in fact, all the functions dehned in (I2.8p are on the same footing, at least for the models 
considered in Section |3l and can only be discriminated by some experiments. 

We want to stress once again that these features were completely hidden in our previous 
analysis, [12], and in many of the papers existing in the literature, since are intrinsically related 
to the explicit appearance of inhnite dimensional Hilbert spaces. 
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